Searching for Supersymmetric Dark Matter- 
The Directional Rate for Caustic Rings. 
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Abstract. The detection of the theoretically expected dark matter is central to par- 
ticle physics and cosmology. Current fashionable supersymmetric models provide a 
natural dark matter candidate which is the lightest supersymmetric particle (LSP). 
The theoretically obtained event rates are usually very low or even undetectable. So 
the experimentalists would like to exploit special signatures like the directional rates 
and the modulation effect. In the present paper we study these signatures focusing on 
a specific class of non-isothermal models involving flows of caustic rings. 

1 Introduction 

In recent years the consideration of exotic dark matter has become necessary in 
order to close the Universe Q. Recent data from the High-z Supernova Search 
Team and the Supernova Cosmology Project || > Q suggest the presence of 
a cosmological connstanta A. In fact the situation can be adequately described 
by a barionic component Qb = 0.1 along with flcDM — 0.3 and — 0.6 (see 
also Turner, these proceedings). 

Since this particle is expected to be very massive, m x > SOGeV, and ex- 
tremely non relativistic with average kinetic energy T < lOOKeV, it can be 
directly detected mainly via the recoiling nucleus. 

Using an effective supersymmetric Lagrangian at the quark level, see e.g . 
Jungman et al jjj and references therein , a quark model for the nucleon fl|Iof 
and nuclear wave functions || one can obtain the needed detection rates. They 
are typically very low. So experimentally one would like to exploit the modulation 
of the event rates due to the earth's revolution around the sun. In our previous 
work [§J ' H we found enhanced modulation, if one uses appropriate asymmetric 
velocity distribution. The isolated galaxies are, however, surrounded by cold 
dark matter , which, due to gravity, keeps falling continuously on them from all 
directions 0. It is the purpose of our present paper to exploit the results of 
such a scenario. 

2 The Basic Ingredients for LSP Nucleus Scattering 

The differential cross section can be cast in the form M: 



da(u, v) 



2(ii r bv) 2 



du 
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a p x o = proton cross-section with i = S, spin, V given by: 

Cp x o — cq (is) 2 (scalar) , (the isovector scalar is negligible, i.e. = a^) 

= a 3 (/° + /A) 2 (spin) , ^ = a (/° + # ) 2 (vector) 
where m p is the proton mass, 77 = m x /mj\fA, and /i r is the reduced mass and 

<t = ^-{G F m N ) 2 ~ 0.77 x 1(T 38 cto 2 (7) 

u = q 2 b 2 /2 or Q = Q u, Q = —L-^ (8) 

where b is (the harmonic oscillator) size parameter, q (Q) is the momentum 
(energy) transfer to the nucleus. In the above expressions F(u) is the nuclear 
form factor and F pp i (u) are the spin form factors [|| (p, p are isospin indices) 
The differential non-directional rate can be written as 

OK ) Yfh 

dR = dR non - dlr = da(u, v)\v\ (9) 

where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the 
detector mass 

The directional differential rate [ fl2| in the direction e is given by : 

dR dlr = £i°l-^ v .eH(v.e) ±- da(u, v) (10) 
m x AniN 27r 

where H the Heaviside step function. The factor of l/2ir is introduced, since we 
have chosen to normalize our results to the usual differential rate. 

We will now examine the consequences of the earth's revolution around the 
sun (the effect of its rotation around its axis is expected to be negligible) i.e. the 
modulation effect. 

Following Sikivie we will consider 2 x N caustic rings, (i,n) , i=(+.-) and 
n=l, 2, ...N (N=20 in the model of Sikivie et al), each of which contributes to the 
local density a fraction p n of the the summed density p of each type i = +, — . 
and has velocity y„ = {y n x,yny,ynz) , in units of vq = 220 Km/s, with respect 
to the galactic center. 



Searching for Supersymmetric Dark Matter 3 



We find it convenient to choose the z-axis in the direction of the motion of 
the the sun, the y-axis is normal to the plane of the galaxy and the x-axis is in 
the radial direction. The needed quantities are taken from the work of Sikivie 
(table 1 of last Ref. by the definitions y n = v n /v Q ,y nz = v n(f> /va,y nx = 

Vnr/vQ,y ny = v nz /v$ . This leads to a velocity distribution of the form: 

N 

f{v>) = J2 S(v - v y„) (11) 

71=1 

The velocity of the earth around the sun is given by ||. 

ve — vq + i>i = vq + v\ ( sina x — cosa cosj y + cosa sirvy z ) (12) 

where a is the phase of the earth's orbital motion, a — around second of June. 
In the laboratory frame we have |§ v = v — ve 

3 Event Rates 

Integrating Eq. (0) we obtain for the total non-directional rate 

R = Rt-^^[1~ h(a,Q min )cosa) (13) 

The integration was performed fromu = m m i n to u — u m ax, where 

u m in = — t: — ,u max = min( —t- , max( — ) , n= 1,2, ...,N) (14) 
Qo a a 

Here y esc = with v escape — 625Km/s is the escape velocity from the galaxy. 
Qmin is the energy transfer cutoff imposed by the detector and a — [v2jU r 6i7o] — 1 . 
Also p n — d n /p,p — Yl n =i dn (f° r each flow +,-). In the Sikivie model |ll| 
2p/p(0) = 1.25. R is obtained (5) by neglecting the folding with the LSP velocity 
and the momentum transfer dependence, i.e. by 

R=^-^y/(^)[Ss + E spin + ^E v ] (15) 

and it contains all SUSY parameters except m x The modulation is described 
in terms of the parameter h. The effect of folding with LSP velocity and the 
nuclear form factor is taken into account by t (see table 1) 

There are now experiments under way aiming at measuring directional rates 
, i.e. the case in which the nucleus is observed in a certain direction. The rate 
will depend on the direction of observation, showing a strong correlation with 
the direction of both the sun's and the earth's motion. In the favorable situation 
the rate will merely be suppressed by about a factor of 2tt relative to the non- 
directional rate. This is due to the fact that one does not now integrate over the 
azimuthal angle of the nuclear recoiling momentum. 
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Table 1. The quantities t and h entering the total non- directional rate in the case 
of the target 53 1 127 for various LSP masses and Q m in in KeV. Also shown are the 
quantities r), h) i = u,d and j = x, y, z, c, s, entering the directional rate for no energy 
cutoff. For definitions see text. 
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We need distinguish the following cases: a) e has a component in the sun's 
direction of motion, i.e. < < tt/2, labeled by i=u (up), b) Detection in the 
opposite direction, tt/2 < 9 < n, labeled by i=d (down). Thus we find : 

1. In the first quadrant (azimuthal angle < < n/2). 

Rdir = R -^j^ r * ~ cosa h i) e z- e 

+ (r'y + cosah 2 + -^-{\cosa\ + cosa))\e y .e\ 

h l 

+ (r* — sinah l 3 + -^-{\sina\ — sffla))|ej;.e|] (16) 

2. In the second quadrant (azimuthal angle tt/2 < <f) < tt) 

Kir = R-^j t^M - cosa h\)e z .e 

+ (r y + cosah l 2 (u) + -^-(\cosa\ — cosa))\e y .e\ 

+ (r l x + sinah l 3 + -^{\sina\ + sina))\e x .e\] (17) 

3. In the third quadrant (azimuthal angle tt < <f) < 3n/2). 

Rdir = R^j - cosa h\)e z .e 

hHu) 

+ (r y — cosah 2 (u) H ^ — (|cosa| — cosa))\e y .e\ 

+ (r l x + sinaHl + (\sina\ + sina))\e x .e\] (18) 

4. In the fourth quadrant (azimuthal angle 3tt/2 < <p < 2tt) 

Rdir = R-^j ^[(rl - cosa h\)e z .e 

+ (r y — cosah 2 + -^-(\cosa\ — cosa)) \e y .e\ 
ti 

+ (r' x — sinoth\ + -^-(\sina\ — sina))\e x .e\] (19) 



4 Conclusions 



We have calculated the parameters describing characteristic signatures needed to 
reduce the formidable backgrounds in the direct detection of SUSY dark matter, 
such as : a) The modulation effect, correlating the rates with the motion of the 
Earth and b) The directional rates, correlated with both with the velocity of the 
sun and that of the Earth (see table 1). 
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We have focused on the LSP density and velocity spectrum obtained from a 
recently proposed non-isothermal model, involving caustic rings pj . Our results 
for isothermal models have appeared elsewhere MM. 

The quantities t and h are given in table 1. We see that the maximum in 
this model does not occur around June 2nd, but about six months later. The 
difference between the maximum and the minimum is about 4%, i.e. smaller 
than that predicted by the asymmetric isothermal models ^,|) . 

For the directional experiments we found that the biggest rates are obtained 
close to the direction of the sun's motion. They are suppressed compared to 
the usual non-directional rates by the factor f rec i — k/(2tt), k = u\. We find 
K ~ 0.7, i — up (observation in the sun's direction of motion) while k ~ 0.3, 
i = down ( in the opposite direction). The modulation is a bit larger than in 
the non-directional case. The largest difference between the maximum and the 
minimum, 8%, occurs not the sun's direction of motion, but in the x-direction 
(galactocentric direction). 

In the case of the isothermal models the reduction factor along the sun's 
direction of motion is now given f re d = to/(Aii t) — n/(2n). Using the values 
of to obtained previously ||, we find that k is around 0.6 for the symmetric 
case and around 0.7 for maximum asymmetry (A = 1.0). The modulation of the 
directional rate depends on the direction of observation. It is generally larger and 
increases with the asymmetry parameter A. For Q m in = it can reach values up 
to 23%. Values up to 35% are possible for large Q m in, but at the expense of the 
total number of counts || . 

Finally in all cases t deviates from unity for large reduced mass. Thus when 
extracting the LSP-nucleon cross section from the data one must divide by t. 
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